Environment-induced sudden transition in quantum discord dynamics 



R. Auccaise/ L. C. Celeri,^ D. O. Soares-Pinto,'^ E. R. deAzevedo,'^ J. Maziero,^ A. 

M. Souza,^ T. J. Bonagamba,^ R. S. Sarthour,^ I. S. Oliveira,^ and R. M. Serra^ 

^ Empresa Brasiletra de Pesquisa Agropecudria, Rua Jardim Botdnico 1024, 

2246O-OOO Rio de Janeiro, Rio de Janeiro, Brazil 
^Centra de Ciencias Naturais e Humanas, Universidade Federal do ABC, 
R. Santa Adelia 166, 09210-170 Santo Andre, Sao Paulo, Brazil 
^Institute de Fisica de Sao Carlos, Universidade de Sao Paulo, 
Caixa Postal 369, 13560-970 Sao Carlos, Sdo Paulo, Brazil 
'^Fakultdt Physik, Technische Universitdt Dortmund, 44^'^! Dortmund, Germany 
''Centra Brasileiro de Pesquisas Fisicas, Rua Dr. Xavier Sigaud 150, 
22290-180 Rio de Janeiro, Rio de Janeiro, Brazil 

Non-classical correlations play a crucial role in the development of quantum information science. 
The recent discovery that non-classical correlations can be present even in separable (unentangled) 
states has broadened this scenario. This generalized quantum correlation has been increasing rel- 
evance in several fields, among them quantum communication, quantum computation, quantum 
phase transitions, and biological systems. We demonstrate here the occurrence of the sudden- 
change phenomenon and immunity against some sources of noise for the quantum discord and its 
classical counterpart, in a room temperature nuclear magnetic resonance setup. The experiment 
is performed in a decohering environment causing loss of phase relations among the energy eigen- 
states and exchange of energy between system and environment, resulting in relaxation to a Gibbs 
ensemble. 



PACS numbers: 03.65.Yz, 03.65.Ud, 03.65.Ta 

The quantum mechanical superposition principle, 
when applied to composite systems, foresees the ap- 
pearance of correlations that cannot be explained in 
a classical context [1]. Initially, the quantum charac- 
ter of a correlated system was attributed to the nonlo- 
cal aspect of quantum mechanics and further associated 
with the violation of Bell's inequalities. The discovery 
of non-separable (entangled) states that do not violate 
the Bell's theorem led eventually to the identification 
of the quantumness of correlations with the separabil- 
ity problem; non-classicality was attributed to entangle- 
ment. Through the development of quantum information 
science (QIS), an operational characterization of entan- 
glement was introduced, as those correlations that cannot 
be generated by local operations and classical communi- 
cation [2]. The development of these ideas led to the 
so-called theory of entanglement, which turned out to be 
a fruitful branch of research [3]. Entanglement is recog- 
nized as an important resource for several tasks in QIS 
[5]. Besides this strong correlation exhibited by entan- 
gled states, there is another kind of non-classical correla- 
tion. A composite quantum system in a mixed state may 
exhibit some non-classical nature in its correlations even 
if it is separable (unentangled) [S]- Such quantum cor- 
relation can be measured as a "gap" between the quan- 
tum and classical versions of mutual information, which 
is the information-theoretic measure of the total correla- 
tion contained in a bipartite system [5]. 

Several approaches have been proposed to quantify this 
generalized quantum correlation present in a bipartite 
mixed quantum state [7| [H]j one of which, quantum 



discord |5j , has received special attention |5HT5] . Beyond 
the importance of non-classical correlations (other than 
entanglement) for the foundations of quantum mechan- 
ics, the relevance of such an issue is increasing in several 
fields as well in applications in QIS. Concerning biologi- 
cal systems, it has been suggested that such correlations 
could play an important role in photosynthesis [9 . In 
condensed matter physics, quantum correlations of sep- 
arable states characterize, even at finite temperature, a 
quantum phase transition |15j . In the context of quan- 
tum field theory, it has been shown that the Unruh effect 
may also lead to an abrupt change in the behavior of cor- 
relations [16] . For applications in QIS, it is interesting to 
know that states with nonzero quantum discord cannot 
be locally broadcast [T7] , and this kind of non-classicality 
is also related to a condition for a complete positive map 
[TT] . Quantum discord was proposed as a figure of merit 
for the quantum advantage in some computational mod- 
els without or with little entanglement [12 [12] • It also 
has a relevant role in mixed-state quantum metrology 

In particular, the quantum discord has been predicted 
to show peculiar dynamics under decoherence [6^. Con- 
sidering two non-interacting qubits (quantum bits) un- 
der the action of Pauli maps, it was shown that, under 
certain conditions, the decay rate of the correlations suf- 
fers an abrupt modification, a phenomenon denominated 
sudden change 0. Moreover, either the classical or the 
quantum part of correlation may be unaffected by deco- 
herence [H [20], giving rise to two distinct decoherence 
regimes, the classical and the quantum. In the present 
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experiment we demonstrate that the aforementioned phe- 
nomena [6l [20] are still present in a real thermal environ- 
ment at room temperature, indicating that such peculiar 
behavior is quite general. We performed such experi- 
mental demonstration in a nuclear magnetic resonance 
(NMR) setup. 

NMR systems have been extensively used as test 
benches for QIS ideas [ITj. Most of these experiments 
were performed in scenarios where the existence of en- 
tanglement was ruled out. The quantum nature of NMR 
systems at room temperature may be ascribed to quan- 
tum correlations of separable states [22]. The main fea- 
ture of the technique for QIS is the excellent control 
of unitary transformations over the qubit provided by 
the use of radio- frequency pulses, which result in unique 
methods for quantum state generation and manipula- 
tion |21| . In general, for room temperature liquid state 
NMR, the Markovian environment in which the qubit 
is immersed can be described by two decoherence chan- 
nels, the amplitude-damping and the phase-damping, 
acting locally on each qubit. Our experiment is per- 
formed on a liquid state carbon-13 enriched chloroform 
sample at room temperature, this sample exhibits two 
qubits, encoded in the and ^'^C spin-1/2 nuclei. The 
two-qubit state is represented by a density matrix in 
the high temperature expansion, which takes the form 
Pab = + eApAB, where e = huiL/AksT ^ 10^^ is 

the ratio between the magnetic and thermal energies and 
ApAB the deviation matrix [3TJ [5S] . 

Several quantifiers of non-classical correlation have 
been proposed. From a thermodynamical approach, Op- 
penheim and coworkers have suggested that the non- 
classical correlations contained in a bipartite quantum 
state can be quantified as the quantum deficit, a differ- 
ence between the amount of work that can be extracted 
from a heat bath, by both and by one of the parts of 
the composite state |8]. In an attempt to build a uni- 
fied framework for correlations, Modi and coworkers have 
proposed a set of measures for the quantification of their 
quantum and classical parts, based on the definition of an 
entropic distance 7 . For the NMR context, one suitable 
generalized non-classical correlation may be computed 
from the experimentally accessible variables as a sym- 
metric version of quantum discord |221 124] (expanded in 
the leading order in e), given by |22| 

QIpab) Hpab) - max I{xab), (1) 
{nj*,nf} 

where the quantum mutual information is given by 

Hpab) « ^ {2tTiApls) - tr [{ApAf] - tr [(Ap^)'] } , 

and the measurement-induced mu- 

tual information is X{xab) ~ 

{2tr [{Axab)'] - tr [{Axa)'] - tr [{Axb?] } , 



with xab = -I- sAxab as the state obtained 

from PAB through a complete projective measure- 
ment map {AxAB = Y.r,j ® nf (ApAB)nf nf ). 

^Pa(b) — ^t^b{a){^Pab} is the reduced deviation matrix 
while Axa(b) stands for the reduced measured deviation 
matrix in the subspace A{B). The classical counterpart 
of Eq. (jlj) is C{pab) max|jj^_jjB|I(xA_B). It is 

worth mentioning that Q{pab) = if and only if pab 
can be cast in terms of local orthogonal basis. In 
other words, the symmetric quantum discord is zero if 
and only if pab has only classical correlations or no 
correlations at all [24]. The aforementioned symmetric 
correlation quantifiers can be computed directly from 
the experimentally reconstructed deviation matrix. 

Let us consider the class of states with maximally 
mixed marginals {pa(b) = Iyi(s)/2), also known as Bell- 
diagonal states: 

PAS - M Iab + ^'"^ ®<^^\^ (2) 

\ i=x,y,z j 

where of^^^ is the ith component of the standard 
Pauli operator acting on the A (B) subspace, a — 
tr [pABf^f (S> <^f] , and Iab is the identity operator. It 
was theoretically predicted [3] that, depending on the ge- 
ometry (encoded in the relations between the parameters 
Ci) of the state ([2]), the evolution of the state correlations 
under the action of Pauli channels presents a peculiar be- 
havior Considering the phase-damping channel, if the 
state components are related as follows: \cz\ > \crc\ , \cy\, 
it was shown the classical correlation is not affected by 
decoherence, while the quantum correlation decays ex- 
ponentially. On the other hand, if \cx\ > \cy\,\cz\ or 
|cy| > |ca;| , Ic^l, and \cz\ ^ 0, the correlations exhibit the 
so-called sudden change in behavior [5]. The classical 
correlation decays exponentially until a specific moment 
in time, thereafter remaining constant, while the quan- 
tum correlation suddenly changes its decay rate at that 
moment. For the other two Pauli channels, we can see 
the same behavior, except that the relations between the 
three components of state ^ are exchanged [5] . 

In order to demonstrate such non-trivial dynamics we 
prepare two initial Bell-diagonal states with specific rela- 
tions between its components. This involves the experi- 
mental mapping of the second term of the right hand side 
of Eq. ([2| onto ApAB- Starting from the thermal equi- 
librium state, the mapping consists of three main steps, 
in which the following are applied to the sample: (i) 
a strongly modulated pulse, (m) a magnetic field gradi- 
ent and (Hi) a pseudo Einstein-Podolsky-Rosen gate im- 
plemented by radio- frequency pulses [21] [25] . After the 
preparation of the state, it is left to evolve freely under 
decoherence and the dynamics of the system is followed 
by quantum state tomography [3] . Figure [l] illustrates 
the pulse sequence in the experimental procedure j27j . 
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FIG. 1. Sketch of the pulse sequence used experimentally 
to follow the dynamics of quantum and classical correlations 
under decoherence. The sequence consists of three blocks: 
the initial state preparation, relaxation delay and read-out 
by quantum state tomography. 



The relaxation process that causes phase decoherence 
and energy dissipation is due, mainly, to internal molec- 
ular or atomic motions that cause random fluctuations 
in the electromagnetic field in which the qubits are im- 
mersed. These fluctuations are characterized by the spec- 
tral densities that encode features of the motion, such as 
geometry and correlation times. The environment can 
be modeled by two independent quantum channels, de- 
scribed as the generalized amplitude-damping and the 
phase-damping channels. In previous investigations of 
sudden-change dynamics [51 [201 US] ) the energy exchange 
channel was not take into account, however in the present 
experiment this relaxation source is unavoidable. It is re- 
markable that the sudden-change dynamics is still hap- 
pen even in the presence of a thermal environment as 
can be seen in Fig. [2| The small decay rate of the classi- 
cal correlation, after the transition point (where it could 
be constant under a phase damping), is due to the pres- 
ence of the thermal noise source. The theoretical pre- 
dictions presented in Fig. [2| are in good agreement with 
the experimental data [57] • We observe a sudden tran- 
sition between two decoherence regimes [B], i.e., classical 
and quantum decoherence [221 . During the first regime 
{t ^ 0.027 s) the decoherence affects more strongly the 
classical aspects of correlation, leaving the quantum as- 
pects with a small decay. After the sudden-change point 
(t ^ 0.027 s) the classical correlation becomes more ro- 
bust against decoherence and the noise degrades more 
effectively the quantumness of correlations. This con- 
firms that the phenomena predicted in [6l [20] for phase 
environments takes place even in the presence of an ad- 
ditional thermal noise. 

Figure [s] also confirms the theoretical prediction [B] 
that, in some states (depending on their geometry), the 
classical correlation may be robust against phase noise. 
The small decay rate of this correlation plotted in Fig. [3] 
comes entirely from the thermal source. Once more, the 
theoretical curves presents very good agreement with the 
experimental data. 

In both of the Figs. 2(c) and 3(c), the theoretical data 
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FIG. 2. (Color online) Sudden change in behavior of correla- 
tions, (a) bar representation of the real and, (b) imaginary 
parts of the initial deviation matrix for the sudden-change ex- 
periment, reconstructed by quantum state tomography. We 
adopted the usual computational basis, where |0) and |1) rep- 
resents the eigenstates of cr^ for each qubit. (c) displays the 
predicted sudden change in behavior of the correlations dur- 
ing their dynamic evolution to thermal equilibrium. The blue 
circles are the experimental data for the quantum mutual in- 
formation, while the red squares and green triangles repre- 
sent the classical and quantum correlations, respectively. The 
black lines are the theoretical predictions. The initial state is 
analogous to the state in Eq. Q with | Ca; | , | Cy | > | | . The 
correlations are displayed in units of (e^/ln2)bit. 



analysis guarantees that the small decay rate presented 
by the classical correlation, where it should vanish, comes 
exclusively from the action of the thermal channel on 
the system. It is worth mentioning that the oscillations 
shown by the experimental curves in both figures are due 
to experimental details (discussed in [371). 

Correlations are ubiquitous in nature. The discovery 
that separable quantum states can exhibit non-classical 
correlation other than entanglement has led to a new un- 
derstanding of the quantum aspects of a physical system. 
Despite the relatively great number of theoretical arti- 
cles concerning decoherence effects, up to date, only few 
experimental investigations of the correlation dynamics 
have been reported previously in literature [551 [23 [H] ■ 
The predictions of refs. [6] [20] were tested in a simu- 
lated phase noise environment |23j . This experiment was 
performed in an optical setup, in which the action of a 
phase-damping channel was simulated in a controllable 
way by a birefringent medium and the qubits were en- 
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UFABC, CNPq, CAPES, FAFESF, and FAFERJ. This 
work was performed as part of the Brazihan National 
Institute of Science and Technology for Quantum Infor- 
mation (INCT-IQ). 
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FIG. 3. (Color online) Immunity against decoherence. (a) 
bar representation of the real and, (b) imaginary parts of 
the initial deviation matrix, reconstructed by quantum state 
tomography, (c) shows that the classical correlation is not 
affected by the action of the phase-damping channel. The 
small decay rate is entirely due to the amplitude-damping 
channel. The displayed pattern is the same as in Fig. 2. For 
this experiment, the initial state is analogous to state in Eq. 
([2| with I I > I Ca; I , I Cj, I . The correlations are computed in 
units of (e^/ln2)bit. 



coded in photon polarizations [23]. Under the action of 
a real global environment, the dynamics of correlations 
was investigated in a NMR quadrupolar system j^ . 

In the present experiment we observed the 
environment-induced sudden change phenomena in 
a real (uncontrollable) thermal noise environment at 
room temperature. The environment-induced sudden 
change takes place in the course of the relaxation of two 
nuclear spins to the Gibbs state. This demonstrate that 
this phenomena may occur in a general context when a 
non-equilibrium system relax to a thermal equilibrium 
state. The methods employed in our experiment to 
follow the dynamics of quantum discord and its classical 
counterpart may be applied to other molecules (in- 
cluding biological ones) and also to other experimental 
contexts where high mixed states are present. The 
two decoherence regimes observed may have important 
consequences for applications in QIS, since the nature of 
correlations in a given mixed state system is somehow 
related to the quantum advantage for performing a 
given task (as for example in quantum metrology |18] ) 
or preventing local broadcasting of information |17j . 
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Supplementary information - Details on 
experimental and calculation procedures 

NMR experiments 

NMR experiments were performed on a two-quantum 
bit system composed of nuclear spins of and ^"^C 
atoms in carbon-13 enriched chloroform (CHCI3). The 
sample was prepared by mixing 100 mg of 99 % ^'Re- 
labelled CHCI3 in 0.2 mL of 99.8% CDCI3 in a 5 mm 
NMR tube (both compounds provided by the Cambridge 
Isotope Laboratories Inc.). NMR experiments were car- 
ried out at 25° C in a Varian 500 MHz Infinity Plus spec- 
trometer located at the Brazilian Center for Physics Re- 
search (CBPF, Rio de Janeiro). A Varian 5 mm dou- 
ble resonance probehead equipped with a magnetic field 



gradient coil was used. 7r/2 pulses of length 7.4 /is and 
8.0 lis were applied to and ^'^C, respectively. For 
specific state preparation, a five-steps strongly modu- 
lated pulses (SMP) H] was designed using a MATLAB 
self-written routine based on a SIMPLEX optimization 
protocol. Spin-lattice relaxation times (Ti) for and 
^^C nuclei, measured by the inversion recovery pulse se- 
quence, were 2.5 s and 7 s, respectively. Spin-spin relax- 
ation times (T2), measured by a CPMG pulse sequence, 
were estimated to be about 1.8 s and 0.29 s for and 
^^C nuclei, respectively. The recycle delay (di) was set 
at 40 s in all experiments. As no refocusing pulse was 
used, the effective transversal relaxation times are given 
by t; = 0.31 s and = 0.12 s for and nuclei, 
respectively. 

The nuclear spin Hamiltonian accounting for the rele- 
vant NMR interactions is given by 



+a;f (if cos if" + If sin tp")+u}f' (if cos ip^ + I^ sin ip^) , (3) 



where if (l^) ^^"^ spin angular momentum operator 

in the a,/? = x,y,z direction for (^'^C); p^ i^p'^^ 
defines the direction of the RF field (pulse phase) and 
uji (wp) is the RF nutation frequency (RF power) for 
(^^C) nuclei. The first two terms describe the Zee- 
man interaction between the "'^H and ^^C nuclear spins 
and the main magnetic field Bg . The corresponding fre- 
quencies are a;///27r 500 MHz and ajc'/27r sa 125 MHz. 
The third term is due to a scalar spin-spin coupling of 
J K, 215.1 Hz. The fourth and fifth terms represent the 
radio- frequency field (RF) that may be applied to the ^H 
and ^'^C nuclear spins, respectively. Besides the interac- 
tions written above, there is a time-dependent coupling 
of the nuclear spins with the environment that includes 
all fiuctuating NMR interactions, which accounts for the 
spin relaxation, e.g., ^H - ^^C dipolar spin-spin couplings, 
interactions with the chlorine nuclei, etc. 

To follow the decoherence of the correlations and dis- 
tinguish the distinct regimes described in the main text, 
an initial state such as that shown in Eq. (1) of the the 
main text needs to be prepared with specific Ci 's . There- 
fore, from the experimental perspective, it is essential to 
prepare a state in which populations and coherences of 
the deviation matrix can be set accordingly. To that 
end, a diagonal state was created from the thermal equi- 
librium by a 4 — 8 ms SMP followed by a 2 ms magnetic 
field gradient pulse. This produced a deviation matrix as 
shown in Eq. Q, where the elements a, /3, 7 and 5 are 
defined by the optimized SMP pulse. Next, this diagonal 
deviation matrix is transformed into an X-type matrix, 
such as that in Eq. ([5]), by a pseudo-EPR gate imple- 



mented by the RF pulse sequence described in reference 



a 
0^00 
7 
5^ 

a + 7 —a + 7 

p+5 -P+5 
+ 5 13 + 6 

—a + 7 a + 7 



(4) 



(5) 



The pulse sequence used in the NMR experiments is 
illustrated in Fig. 1 of the main text. The SMP was 
designed to produce a diagonal state with populations 
a, /3, 7 and 5, so that, after the transformation by the 
pseudo-EPR gate, a state with a deviation matrix in the 
form equivalent to the state in Eq. (2) of the main text, 
with specified Ci's, was obtained. After this state prepa- 
ration step, the system was left to evolve for a period 
T = m/4J, with m = 1,2,3,4, ... During this evolution 
period (considering on-resonance evolution), two terms 
of the Hamiltonian shown in Eq. ^ act on the nuclear 
spins: the scalar coupling, which produces oscillations of 
the quantum coherence, and the time-dependent term, 
which accounts for the system relaxation. It is worth 
mentioning that the elements in the main and cross di- 
agonals of the deviation matrix account for populations 
and quantum coherences that do not oscillate with the J 
coupling evolution, so that, for a perfect X-type devia- 
tion matrix, the aforementioned oscillation should not be 
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observed. However, the experimental deviation matrices 
have small coherences off the X positions, which oscillate 
during the evolution. This leads to small oscillations in 
the spectral line intensities used to perform the state to- 
mography and, because the tomography process involves 
solving coupled equations, this oscillatory behaviour is 
captured in the main and cross diagonal elements, caus- 
ing the oscillations observed in the correlations shown in 
Figs. 2 (c) and 3 (c) of main text. The oscillatory evo- 
lution may be suppressed from the experimental data by 
setting the evolution step to a multiple of 1/ J . However, 
in order to show unequivocally the sudden change in the 
quantum correlations, the evolution step was set to 1/4 J, 
giving rise to the small oscillations observed in the evolu- 
tion of the correlations depicted in Figs. 2 (c) and 3 (c) of 
the main text. After the evolution period, the quantum 
state was read-out by quantum state tomography, as de- 
scribed in reference O |4] . Multi-step increments of r in 
successive experiments allowed the effect of the spin en- 
vironment on the initial deviation matrix to be followed. 
Final states corresponding to 251 distinct r values were 
acquired. The high homogeneity of the static magnetic 
field was guaranteed by the linewidth of w 0.95 Hz in ^H 
spectra. 



Decoherence analysis 

The decoherence process is theoretically analysed 
through the operator sum representation technique 
in which the evolution of the density operator is given by 

p(t) = ^ii;,(t)p(o)4 (t), 

k 

where Ek {t) are the well-known Kraus operators. 

Amplitude- Damping Channel: In this scenario, the 
amplitude-damping channel is described by the follow- 
ing set of operators 



where, in the NMR context, 7 = 1/2 — e/2 and p = 
1 — exp (— t/Ti), Ti being the longitudinal relaxation time 
of the qubit under consideration. We observe that in our 
case, the relaxation times are different for the two qubits, 
since they have distinct Larmor frequencies. 

Phase-Damping Channel: For this case we have 

where A = 1 — exp (—t/T2) with T2 being the transverse 
relaxation time associated with the qubit. 

The agreement between the model and the experimen- 
tal data, as observed in Figs. 2 (c) and 3 (c) of the 
main text, is fairly good. In the theoretical model we ne- 
glected a number of possible error sources, such as mag- 
netic field inhomogeneities or oscillations of the residual 
off-diagonal elements of the deviation matrix. In spite of 
that, all the singular features predicted by the quantum 
model are very clear in the experimental data. 

Within this model and taking into account that both 
environments act independently on each qubit, it is 
straightforward to eliminate from the experimental data 
the contribution of the amplitude-damping channel to 
the dynamics of the correlations. This procedure removes 
the small decay rate of the classical correlation (which it 
should be constant) in both Figs. 2 (c) and 3 (c) of the 
main text, proving that this correlation is really immune 
to the phase-damping channel in certain cases, corrobo- 
rating the theoretical predictions made in ref. |6]. 
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